AH  6l>1)  4  06 


WAVE  FORCES  ON  PILES:  A  I) 1 1  E  K  A  C  I  I  O N  I'll  I-  OK  Y 
R .  (’ .  M  a i’Ca  m  y ,  cl  a  I 

Corps  of  Engine o  i  s 

Washington,  I).  O. 

December  ll>54 


This  document  has  been  approved  for  public  release  and  sale. 


OFFICE  OF  THE  CHIEF  OF  ENGINEER* 


WAVE  FORCES  ON  PILES 
A  DIFFRACTION  THEORY 


TECHNICAL  MEMORANDUM  NO.  M 


Nit  COPY 

BETTIES  TO 
•EACH  EROSION  BOARD 


?“  document  has  been  appro, 
tor  Public  roleose  and  o.Ho; 
H»  distribution  in  unllnltod 


WAVE  FORCES  ON  PILES: 
A  DIFFRACTION  THEORY 


r. 


TECHNICAL  MEMORANDUM  NO.  69 
BEACH  EROSION  BOARD 
CORPS  OF  ENGINEERS 


DECEMBER  1954 


FOREWORD 


Although  circular  piling  is  a  much-used  structural  element  in  shore 
protection,  harbor,  and  other  maritime  structures,  only  recently  have 
significant  advances  been  made  toward  gaining  a  quantitative  understand¬ 
ing  of  the  forces  developed  by  wave  action  against  piling.  The  present 
report  deals  with  this  subject. 


The  report  was  prepared  at  the  University  of  California,  Berkeley, 
California.  The  work  on  which  the  report  is  based  was  sponsored  by  the 
Office  of  Naval  Research,  U.  S.  Department  of  the  Navy.  The  authors  of 
the  report  are  R.  C.  MacCamy  and  R.  A.  Fuchs  of  the  Institute  of 
Engineering  Research,  University  of  California.  Because  of  its  appli¬ 
cability  to  the  research  and  investigation  program  of  the  Beach  Erosion 
Board,  and  through  the  courtesy  of  the  authors,  the  report  is  being 
published  at  this  time  in  the  technical  memorandum  series  of  the  Beach 
Erosion  Board.  Views  and  conclusions  stated  in  the  report  are  not 
necessarily  those  of  the  Beach  Erosion  Board. 
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WAVE  FORCES  ON  PILES:  A  DIFFRACTION  THEORY 


by 

R.  C.  MacCamy  and  R.  A.  Fuchs 


Introduction.  This  report  contains  two  main  results.  In  the  first 
section  an  exact  mathematical  solution  is  presented  for  the  linearized 
problem  of  water  waves  of  small  steepness  incident  on  a  circular  cyl¬ 
inder.  The  fluid  is  assumed  to  be  frictionless  and  the  motion  irro- 
tational.  This  section  includes,  in  addition  to  the  formal  mathematical 
treatment,  some  simple  deductions  based  on  the  assumption  of  very 
small  ratio  of  cylinder  diameter  to  incident  wave-length.  The  prin¬ 
cipal  results  of  the  theory  are  summarized,  for  convenience  in  calcu¬ 
lations,  in  the  second  section.  Also  presented  are  none  suggestions 
as  to  possible  extensions  of  the  theory  to  take  care  of  more  extreme 
wave  conditions  and  other  obstacle  shapes. 

*  -• 

The  second  result  is  an  attempt  to  apply  the  theory  to  the  com¬ 
putation  of  actual  wave  forces  on  cylindrical  piles.  The  basis  of 
comparison  is  a  series  of  tests  performed  in  the  wave  channel.  The  .. 
agreement  is  found  to  be  quite  good  in  the  region  in  which  the  as¬ 
sumptions  of  the  theory  are  fairly  closely  realized. 

Theory.  The  problem  of  diffraction  of  plane  waves  from  a  circular 
cylinder  of  infinite  extent  has  been  solved  both  for  electromagnetic, 
and  sound  waves.  Only  alight  modifications  are  necessary  to  obtain  a 
corresponding  solution  for  water  waves  incident  on  a  circular  pile. 
Reference  is  made  to  Morse (l)*  especially  for  the  expansions  in 
equations  2,  3,  $,  and  for  a  survey  of  the  asymptotic  developments  of 
the  Bessel's  Functions. 

The  following  assumptions  are  made.  The  fluid  is  frictionless 
and  moving  irrotationally.  The  ratio  of  the  height  of  the  waves  to 
their  length  is  sufficiently  small  so  that  all  quantities  involving 
the  parameter  (H/L)  in  the  second  or  higher  powers  may  be  neglected 
without  sensible  error,  thus  giving  rise  to  the  so-called  linear 
theory.  The  waves  are  indident  on  a  vertical  circular  cylinder  which 
extends  to  the  bottom.  The  depth  of  the  water  is  d,  finite.  > 

A  set  of  axes  x  ,  y,  z  is  chosen  with  z  directed  positively  up¬ 
ward  from  the  still-water  level.  The  cylinder  of  radius,  a,  is 
assumed  to  lie  along  the  z-axis  and  cylindrical  waves  ate  incident 
from  the  negative  x-direction.  The  velocity  potential  of  the  inci¬ 
dent  wave  then  may  be  written, 

<*>(*)-  2£sh.k  i  (kx-crt)  (l) 

“  Ter  cosh  kd  v  ' 


♦Numbers  in  parentheses  refer  to  list  of  references  on  page  11  . 
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It  is  understood  here  that  the  actual  potential  is  the  real  part  of 
this  oomplex  expression,  and  that  in  order  to  find  the  physical 
solution  in  what  follows,  it  is  necessary  to  te*ce  real  parts* 


Introducing  polar  co-ordinates  r  and  0,  equation  1  admits  of 
an  expansion  in  cylindrical  harmonics,  having  the  formt 


-f  Got*)  -  S**- « •  w 


The  assumption  is  now  made  that  the  reflected  wave  admits  of  a 
similar  expansion*  The  particular  combination  appropriate  to  a  wave 
moving  outward,  symmetrically  with  respect  to  0,  that  is  such  that 

f  (•  ©)  •  *  (©),  i», 

^ArOorhO  {jn(kr)  ♦  i  Tm(krjJ-  +i9t  (3) 

This  combination  of  Bessel  Functions  is  known  as  the  Hankel  function 
of  the  first  kind,  H^'  (kr),  and,  for  large  values  of  r,  has  the 
asymptotic  form 

H^1)  (kr e1^  *  T  ^  (1*) 


Henoe  equation  3  has,  for  large  values  of  r,  the  form  of  a  periodic 
disturbance  moving  outward  in  the  r  direction,  with  frequency  <r  and 
wave  number  k,  and  vanishing  at  r  aCD  . 

For  the  total  velocity  potential,  0,  there  is  taken  a  super¬ 
position  of  0  (i)  end  an  infinite  series  of  terms  like  the  quantities 
Ag  ere  then  determined  by  setting  the  particle  velocity  normal  to  the 
oylinder,  that  ia  ,  equal  to  zero  at  the  surfaoe,  r  ■  a* 


The  result  of  this  calculation  is, 

*  •  U  •lwt  [J°(kr) 


Jo1  Oca) 
H0^)’  (ka) 


H0(2)  (kr) 


<2I,  ■  ^y(k.)H.(8)  (kr))  eo*  "*]  (5) 

where  Hlt(2)(kr)  is  the  Hankel  Function  of  the  second  kind  and  equals 
J-*i  Tn*  This  result  is  given  by  HavelockC?)  for  the  special  case 
of  Infinite  depth. 
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The  pressure  exerted  on  the  cylinder  is  computed  from  Bernoulli's 
equation, 


mP$t  --f  (d  ~)z\-  <’a  ^2-gfZ 
■ F a t  r.  'a  /  '  ay-  y^ 


P  [d  1  a  y'  "yAi  (6) 

where,  in  the  li^ar  theory,  the  squared  terms  are  neglected. 

The  x -component  of  the  force,  per  unit  length  in  the  z-direction, 

rn 

Fz  ■  Re  2  I  p  (9)  a  cos  (tt  -  0)  d9 

-*Ci 


is, 


Only  the  term  in  cos  9  will  contribute  to  this  integral  and  the  result 
after  taking  the  real  part  may  be  written  as, 


(7) 


where 


J,'  (ka) 
tanO  s  - - 

V  (ka) 


A  (ka) 


4V2  (ka)  ♦  Yx'2  (ka) 


These  functions  are  plotted  in  Figures  1  and  2,  ka  being  equal  to 
TT  DA* 

The  moment  about  a  point  z  «  u,  on  a  cylinder  extending  to  depth 
v  below  the  still -water  level  may  be  easily  computed  from  equation  7, 
assuming  that  the  motion  of  the  fluid  is  the  same  as  if  the  cylinder 
extended  to  the  bottom.  The  expression  for  the  moment  is, 


"Hi 


v  -  f  (z  -  u) 

*  •<_  u 


Fz  dz 


(8) 


To  be  cc.isistent  with  the  linear  theory  the  integration  need  only 
be  carried  up  to  the  still-water  level  z  s  0,  the  result  being 

m.  m  2g/>  H  [" uk  sinh  kd-sinh  k(d-v)-vk  sinh  k(d-y)+cosh  kd 

"  "  iT  K  '  L  cosh  kd 

cos  (at  -a)  (9) 


cosh  k(d-v) 
cosnkd 


The  special  case  of  a  pile  hinged  about  the  bottom  is  evaluated  by 
setting  u  s  -d,  v  =  d. 
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^  .  Lisp  *(k.)  (K  d  «V»fc  kd  -  V  con  (<r  t  -a)  (10) 

The  function  D  (kd)  ■  1  -  cowh  kd^»  kd  slnh  kd  giving  the  dependence 

on  depth  is  plotted  in  Figure  3* 


An  estimate  of  the  effect  of  second  order  terms  on  the  moment  m^y 
may  be  immediately  obtained  from  equation  8  by  evaluating  that  portion 
of  the  integral  from  zero  to  f)  .  To  the  second  order,  for, 

jf  m  ^  sin  a  t  (13) 

Xi)  2 

(a  -  u)  Fa  da  •  £iLj|Lj}  A(ka)  sin«rt  cos  (et  -a  )  (12) 


This  calculation  omits  that  portion  of  the  second -order  terms  arising 
from  the  second  term  in  the  velocity,  but  this  latter  term  may  be  ex¬ 
pected  to  be  small.  It  is  noted  that  the  result  (12)  may  be  obtained 
by  assuming  that  the  force  and  lever  arm  are  constant  over  the  range 
having  the  value  at  t  ■  0  and  multiplying  these  constant 
values  toy  the  length,!)  .  For  the  special  case  of  a  cylinder  hinged  at 
the  bottom  the  total  moment  becomes 


& 


A(ka)  D(kd)  cos^t  1  ♦ 


1 


k*Hd  Bln. 


■] 


(13) 


From  equation  13  it  is  seen  that  the  maximum  moment  occurs  for, 

|2 


sin  (r  t^,v 


**  LhrarJ 


(1U) 


k  Hd 


c  ihtot 

and  has  the  value  obtained  by  substituting  (#*t)  into  equation  13. 


For  cylinders,  the  diameters  of  which  are  small  compared  to  the 
length  of  the  waves,  the  foregoing  theory  admits  of  several  simplifi 
cations.  Asymptotic  values  of  the  Bessel's  Functions  and  their 
derivatives  are  presented  for  reference  in  Table  IV.  These  lead 
immediately  to  the  approximate  formulas, 


A(ka)  (ka)2 

a(k.)?J-  (k.)2, 


(16) 
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In  particular,  aquation  7  m/  than  ba  ra  pi  a  cad  by 


.  .2  cosh  k  (d*e) 

f»  m*  r  *  "  k  *  co.l'  ii  ‘ 


coa  r  t 


In  thin  fora  tha  forca  Fg  adaita  of  a  auch  simpler  derivation. 

For  a  van  incident  on  a  vortical  wall  at  an  arbitrary  angle  there  ia 
complete  reflection  without  loaa  of  energy,  resulting  in  a  total 
praaaura  equal  to  twice  that  of  tha  incident  wave.  Assuming  that  thia 
reault  hold*  for  tha  cylinder  alao,  an  incident  wave  with  velocity 
potential  given  by  equation  1  will  give  riae  to  a  real  pro  Mure , 

pe-fgH  coah  h  (dea)  ain(kx-^t).  (17) 

co ah  Yd 

The  resulting  forca,  F,  la  than  obtained  by  integration  aa  for 
equation  7,  giving  the  relatlonehip 


•2  g  H 


coah  k  (deal 


J  aln  (ka  coa  9  -  - 1)  coa  9  d9 


But  no  Tor  Mall  valuea  of  ka,  expanding  the  integrand  in  equation  18 
given 


v  ~ 
Fa  ■ 


U  »f  g  H 


coah  k  (d*aj 


coah  t  t  ka  coa*9  d9,  (19) 


which  leada  again  to  equation  7*.  It  ia  to  be  noted  in  connection 
with  thia  equation,  that  the  force  1^  ia  eaual  to  the  ao -cal led 
"virtual  naaa  force"  in  Morlson's  reault**/  provided  the  experimental - 
ly  determined  eonataat  Cm  ia  taken  aa  two.  The  result  la  to  be  ex¬ 
pected  since  an  essential  assumption  of  Horiaon'a  theory  la  that  the 
form  of  the  incident  wave  ia  little  affected  by  the  presence  of  the 
cylinder.  From  equation  2}  it  la  seen  that  thia  aaaumption  ia 
equivalent  to  the  smallness  of  the  ratio  of  pile  diameter  to  wave 
length.  It  ia  to  be  noted  in  this  connection  that  the  exact  theory  of 
the  present  report  repreaente  an  extension  since  its  accuracy  does 
not  depend  on  the  relative  site  of  the  cylinder.  The  value  of  Cg 
quoted  by  Morlaon  for  a  aerlea  of  model  studies  ia  nearly  1.5. 

Thia  type  of  analysis  admits  of  certain  extensions.  For  example 
the  same  technique  might  be  used  to  obtain  foroea  on  more  oompl'cated 
shapes,  the  dimensions  of  which  are  mull  compered  t c  the  wave  «ength, 
since  a  knowledge  of  the  form  of  the  reflected  wave  ia  not  neoessary. 
It  ia  alao  shown  in  the  next  section  how  an  estimate  of  the  effect 
of  steeper  waves  may  be  obtained  in  a  similar  manner. 


A  more  exact  analysis  of  the  relative  effects  of  the  incident 
and  reflected  waves  in  possible  from  the  small  cylinder  theory  and 
will  offer  justification  for  the  developments  of  the  preceding  para¬ 
graphs.  The  surface  profile  may  be  obtained  from  the  velocity 
potential ,  0,  given  by  equation  5  from  the  formula, 

Jim  I  l  .  o  (20) 

n  *  at 


thin  gives 
a  a 


OO 

♦2  I  in 

e«i 


1  I 

n — —  cos  n6  I 

'  Oca)  J 


where  use  has  been  made  of  the  identity, 

Jm(x)  (*>  "  V  <x>  «m(2)  (*)  ■  “  Jrj  (22) 

Using  the  asymptotic  formulas  for  the  Bessel  Functions  for  small 
values  of  ka,  equation  21  becomes,  on  taking  the  real  part, 

(1))r«aV  £  +  U  (ka)^  cos^O  sin  (rt-f)  (23) 


where 


tan  ^  s  2ka  con  9 


In  the  same  notation  the  pressure,  at  the  surface  of  the  pile,  in, 
r  —  ~\  oosh  k  (d*z) 


&4ii  - 1 -  *2  ?  in - - -  cos  n  9  cofl> 

Vks  /2V  Ml  (2) « 

_HqU;  (ka)  Hn'  '  (ka)  J 


irka  /2)'  Ml  (2)« 

[V  (k.)  '  (ka)  J  (21*) 

or  for  small  piles,  the  real  part  of  equation  21*  gives 

P  "s  gPH  (ninO’t  ♦  2ka  cos  0  costf’t)  °°-8^ (25) 

cosh  ltd 

It  can  easily  be  shown  that  the  pressure  due  to  the  incident  wave 
only  is  to  the  same  degree  of  approximation, 

a  gP H  (sinat  ♦  k  a  coa  0  cos  r  t)  ^cosK'  ^  C2^) 

so  that  the  pressure  due  to  the  reflected  wave  is, 

P(l)  •  .OH  ■»  CO,  .  .08. t  *£* W 
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It  is  observed  that  the  first,  and  largest,  terra  of  equation  26 
is  independent  of  6  and  hence  will  contribute  nothing  to  the  force, 
Fz.  Hence  the  "effective"  pressures  due  to  the  incident  wave  and 
the  reflected  waves  are  identical.  This  is  in  contrast  to  the  effect 
on  the  surface  elevation,  since  equation  23  shows  that  the  deviation 
from  that  of  the  incident  wave  alone  is  small. 


Summary.  The  diffraction  of  long-crested  waves  incident  on  vertical 
circular  cylinders  extending  from  above  the  water  surface  to  the 
bottom  is  treated  exactly  within  the  framework  of  the  linearized 
irrotational  theory.  The  essential  results  are  summarized  below. 

Letting  z  be  the  distance  along  the  cylinder,  in  the  direction 
of  its  axis,  with  positive  direction  upward  from  the  still-water 
level,  the  x-component  of  the  force  on  the  cylinder  per  unit  length 
in  the  z-direction  and  at  depth  z,  is, 

Tz  =  l£ji-S  A(g)  CO.  (<Tt -a  )  (28) 

Jl'  (ir£  ) 

ri'  (,rr  > 

A  (r)  ■  .  ■  ■■■■■■  i-  I~ 

'Jv  ('?)  »v  (wt> 

when  the  surface  elevation  is  given  by, 

T)*  ^  sin  (k  x  -  9  t)  (29) 

and  Y^  are  t*1®  Bessel's  Functions  of  the  first  and  second  kinds, 
respectively,  and  primes  indicate  differentation.  The  functions  A 
and  (1  are  plotted  in  Figures  1  and  2.  Additional  values  can  be  ob¬ 
tained  from  a  set  of  tables  published  by  the  Mathematical  Tables 
Project(b). 


where 

tan  OL  • 


The  corresponding  movement  on  a  cylinder  extending  to  depth  v 
below  the  still -water  level  and  hinged  at  depth  u  relative  to  the 
still-water  level  is  given  by 

-2b PH  . ,Dx  I  uk  sinh  kd  -  sinh  k(d-v)  -  vk  sinh  k(d-v) 
mu'v  «  -jf-  A(r)  - — 


♦ 


co^h  kd  -  cosh  k(d-v) 
cosh  kd 


cos  (<T  t  -a  ) 


(30) 
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In  the  special  case  of  a  cylinder  extending  to  the  bottom  and 
hinged  at  the  bottom,  u  ■  -  d  and  v  =  d  and  equation  30  becomes 


where 

D  (kd)  ■  1  ~  cosh  kd  +  kd  sinh  kd 

"  cosh  kd 

The  function  D  (kd)  is  plotted  in  Figure  3*  The  moments  in  this  case 
may  be  easily  computed  through  the  use  of  Figures  U  and  5.  Assuming 
H#  T,  d,  D  to  be  known,  the  ratio  D  is  found  from  Figure  U  and  then  mo 
computed  from  Figure  5.  I 


For  the  case  of  small  cylinders,  that  is,  such  that  the  ratio 
of  the  diameter  to  the  wave  length  is  small,  these  formulas  may  be 
greatly  simplified.  This  appears  to  be  the  most  important  case  as 
is  seen  by  considering  Figure  2.  For  a  l£0  foot  ocean  wave,  the 
cylinder  diameter  could  exceed  fifteen  feet  without  appreciable  de¬ 
viation  from  the  approximate  formulas.  For  this  condition  the  functions 

A  (^)  ondd(~)  may  be  replaced  by 


A(S) 


-5-  ) 


-  /D  \  <v  T* 

a  lr  '  =  T" 


<E>2 


The  force,  Fz,  then  becomes 

F_  -  ,T 2  £.  K  -p2  (  H  >  cosh  k  (d  ♦  z)  coso-t  (28') 

z  2  L  cosh  kd 

and  the  surface  elevation  at  the  circumference  of  the  pile  may  be 
written, 

*)  z  j  \j  1  ♦  cos2  9  sin  (<Tt  -  f  )  (33) 

where 

ten  yr  coa  9  , 


while  the  pressure,  at  depth  z,  as  a  function  of  0,  is 


(3U) 


8 


where 


cos  9 


tan 


2ffD 

L 


A  comparison  of  equations  29  and  30  indicates  the  maximum  force  and 
moment  occur  almost  ninety  degrees  out  of  phase  with  the  crest  of  the 
wave,  that  is  approximately  at  the  time  the  wave  is  passing  through 
the  still-water  level. 


In  the  formulas  thus  far  presented  the  linear  theory  has  been 
strictly  followed.  Approximations  to  the  effects  of  steeper  waves  may 
be  obtained  by  making  some  additional  assumptions.  It  has  been  shown 
previously  that  in  the  case  of  small  piles  the  force,  Fz>  given  by 
equation  7'  is  exactly  twice  that  of  the  incident  wave  alone.  Assuming 
that  this  is  a  general  result,  the  second  and  higher  order  terms  in 
the  parameter  (")  may  be  introduced  into  the  force  calculations.  To 

L 

the  second  order,  the  force  obtained  in  this  manner  is, 


2  F. 


P  TTg  D  d*  * 


cosh 
cosh 


sh  kd  1  1  \  sinh-hcd  cosh  kd 


2  sinh  2  kd 


(24£> 


sin  2  9  t  ( 3U ) 


for  the  surface  elevation  of, 

^/H  =  1_  cos  a  t  ♦  7-  v  (ji)  ctnh  kd  (1  *  — 
2  U  L  2 


sinh^.  ’  kd 


-)  cos  2  <r  t 


(35) 


For  purpose  of  calculation  a  set  of  force  distribution  curves 
has  been  presented  in  Figures  6  and  7,  The  corresponding  moments  may 
be  computed  graphically  according  to  the  following  procedure.  For 
the  moment  about  a  hinge  at  depth  compute  Zi/d  on  the  vertical 
scale.  A  new  curve  then  may  be  plotted  with  abscissa  z  z, 

’  r)d 

times  the  old  abscissa,  and  the  corresponding  moment  will  be  equal  to 
the  area  under  this  curve,  after  multiplication  by  the  respective 
numerical  factors.  The  coefficients  of  the  sin  <Tt.  and  sin  2  t 
terms  in  the  force  equation  have  been  designated  F^l)  and  Fz(2), 
respectively. 

The  finite  height  of  the  waves  introduces  a  second  correction  to 
the  calculated  moments,  namely  the  contribution  to  the  total  moment 
of  that  portion  of  the  wave  above  or  below  the  still-water  level. 

For  a  pile  hinged  at  position  u  this  correction  term  is, 
approximately, 
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(36) 


2 

s  — Jl  1! — ii_  A  (  ^  )  cos  (ft  -  a)  sin  <T  t 
u  k  L 

Comparison  with  Experiment.  A  series  of  experiments  has  been  carried 
out  by  Morison(b)  in  the  wave  channel  to  measure  moments  on  cylindrical 
piles  under  varying  sets  of  wave  conditions.  The  cylinders  were 
hinged  at  varying  depths  and  subjected  to  regular  wave  trains  which 
were  of  essentially  three  types;  moderately  steep  waves  in  shallow 
water,  steep  waves  in  deep  water  and  low  waves  in  deep  water.  In  Table  I 
are  presented  the  results  of  these  experiments  for  piles  hinged  on 
the  bottom  in  low  waves  in  deep  water.  The  theoretical  moments,  m0, 
computed  from  the  graphs  of  Figures  U  and  5,  are  corrected  for  the 
finite  height  of  the  waves  by  adding  as  given  by  equation  36. 

In  Table  II  the  results  for  the  largest  cylinder  in  the  same 
wave  conditions  are  presented  w ith  the  pile  hinged  at  varying  depths, 
z.  It  is  seen  that  in  both  of  these  tables,  in  which  the  actual 
conditions  approximate  the  assumptions  made  in  solution,  the  agreement 
is  good. 

In  Table  III  the  results  for  the  first  two  types  of  waves  are 
presented.  The  first  three  entries  correspond  to  moderately  steep 
waves  in  shallow  water,  and  the  last  three  to  steep  waves  in  deep  water. 
The  deviations  here  are  seen  to  be  quite  large,  reflecting  the  fact 
that  the  waves  cannot  be  closely  approximated  by  sine  waves  in  this 
range  of  H  and  d  . 

L  L 

Conclusions  and  He  commendations  for  Further  Work.  The  rather  large 
deviations  of  the  experimental  results  from  calculated  values,  which 
are  indicated  in  Table  III,  give  rise  to  the  need  for  a  consideration 
of  possible  sources  of  error  together  with  possible  modifications. 

In  order  to  obtain  agreement  with  experiment  MorisonO)  has  intro¬ 
duced  a  second  component  of  force  on  the  pile  which  he  designates  as 
a  "drag"  force.  It  has  been  previously  pointed  out  that  hir>  accelera¬ 
tive  force,  in  the  special  case  of  small  piles,  may  be  identified 
with  the  diffraction  theory  of  this  report,  provided  C_  is  taken 
equal  to  two.  The  introduction  of  the  drag  force  is  tnen  equivalent 
to  the  assumption  that  drag  and  diffraction  forces  may  be  separated, 
each  being  considered  to  act  independently  of  the  other,  an  assumption 
which  may  not  be  well  justified. 

The  force  attributed  to  "drag"  is  essentially  of  two  parts.  One 
arises  from  the  viscosity  of  the  fluid  and  the  corresponding 
frictional  drag  exerted  by  the  fluid  moving  past  the  cylinder.  This 
problem  has  been  considered  approximately  using  Schlicting's  theory 
of  periodic  boundary  layers  and  the  results  indicate  that  frictional 
efrects  are  unimportant.  The  second  part  of  the  drag  force  is  due 
to  the  separation  of  the  lines  of  flow,  with  the  resultant  decrease 
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in  pressure  behind  the  cylinder.  The  wake  behind  the  cylinder  is 
then  essentially  a  region  of  no-motion  except  for  the  possible  forma¬ 
tion  of  vortices.  The  exact  nature  of  this  wake  is  not  well  under¬ 
stood  even  for  the  case  of  steady  flow  and  very  little  is  known 
about  periodic  motion,  for  then  there  is  continual  change  with  in¬ 
creasing  and  decreasing  velocity. 

For  the  diffraction  theory  presented  in  this  report  the  motion  is 
symmetrical  around  the  cylinder.  Hence  no  separation  occurs  and  the 
wake  drag  must  be  zero.  In  light  of  this  result  it  seems  doubtful  that 
the  correction  due  to  drag  could  be  made  simply  by  addition  of  a  term 
corresponding  to  a  wake  while  still  maintaining  the  same  diffraction 
force. 


Drag  forces  are  determined  experimentally  for  the  case  of  steady 
flow  past  a  cylinder  in  the  following  manner.  The  assumption  is  made 
that  the  drag  force  is  proportional  to  the  square  of  the  velocity, 
the  diameter  of  the  cylinder  and  the  density  of  the  fluid.  The  con¬ 
stant  of  proportionality,  called  Cp,  is  then  determined  empirically 
for  various  valuer  of  the  Reynolds  numuer.  *'*orison  has  assumed  that 
this  result  will  also  hold  for  periodic  motion,  an  assumption  which 
needs  considerable  investigation,  since  the  flow  behind  the  cylinder 
may  not  be  able  to  adjust  rapidly  enough  to  maintain  steady  state 
conditions.  Some  calculations  have  been  made,  however,  using  Moriscn' s 
assumptions  and  it  is  found  that  the  introduction  of  drag  does  not 
appreciably  improve  the  results. 

The  results  of  this  report  indicate  that  a  great  deal  of  additional 
work  might  profitably  be  carried  out.  In  particular  a  detailed  ex¬ 
perimental  study,  with  photographs,  of  the  actual  state  of  motion  be¬ 
hind  the  cylinder  would  be  of  considerable  value  in  estimating  the 
effect  of  drag.  It  nicht  be  expected  that  for  moderately  snail 
velocities  the  motion  up  to  a  certain  point  on  the  cylinder  is  well 
approximated  by  the  diffraction  heory,  while  beyond  that  point  the 
flow  separates  leaving  a  dead  water  region.  If  this  should  prove  to 
be  the  case,  additional  theoretical  results  are  possible. 
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•  TABLE  I 


D 

in.) 

d 

r 

H 

r 

(n^)  exp. 
(ft.  lba.) 

(uiq)  theo. 
(ft.  lbs.) 

mo  ♦Amo 
(ft.  lbs. 

* 

0.1*0 

0.037 

0.0207 

0.0202 

0.0205 

0.1*1 

0.038 

0.0203 

0.0207 

0.0211 

1 

0.1*0 

0.036 

0.0903 

0.0801* 

0.0816 

i 

0.1*0 

0.037 

0.0998 

0.0813 

0.0825 

2 

0,1*0 

0.037 

0.2910 

0.3U* 

0.320 

2 

0.1*0 

0.037 

0.2905 

0.310 

0.315 

TABLE  II 
D  •  2  inches 


z 

2 

d 

r 

H 

r 

("o)  exp. 
(ft.  lbs.) 

(n.)  theo. 
(ft?  lbs.) 

(ny,  ♦Amo) 
(ft.  lbs.) 

0.25 

0.39 

0.037 

0.0335 

0.0311 

0.0373 

0.1*2 

0.1*0 

0.039 

0.0836 

0.0808 

0.0891* 

0.52 

0.1*0 

0.037 

0.112 

0.110 

0.120 

0.68 

0.1*0 

0.038 

0.151* 

0.176 

0.187 

0.78 

0.39 

0.037 

0.205 

0.221 

0.232 

0.98 

0.1*0 

0.037 

0.291 

0.315 

0.320 

0.98 

0.1*0 

0.037 

0.291 

0.310 

0.315 
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TABLE  III 


D 

(in.) 

d 

r 

H 

L 

(mo)  exp. 
(ft.  lbs.) 

(u^j)  theo. 
(ft.  lbs.) 

(m  ♦  A  m 
(ft.  lbs 

i 

0.15 

0.066 

0.1158  • 

0.0650 

0.0703 

1 

0.16 

0.066 

0.0785 

0.0603 

0.0628 

1 

0.39 

0.093 

0.3900 

0.2065 

0.2260 

2 

0.16 

0.065 

0.2205 

0.1530 

0.1610 

2 

0.60 

0.090 

0.962 

0.822 

0.886 

TABLE  IV 


Asymptotic  Expansions  for  Bessel' s  Functions  and  Their 
Derivatives  for  Small  x. 
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